Total Variation of the Control and Energy of Bilinear Quantum Systems 
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Abstract — In the present note, we give two examples of 
bilinear quantum systems showing good agreement between the 
total variation of the control and the variation of the energy 
of solutions, with bounded or unbounded coupling term. The 
corresponding estimates in terms of the total variation of the 
control appear to be optimal. 

I. INTRODUCTION 

A. Control of quantum systems 

The state of a quantum system evolving in a Riemannian 
manifold Q, is described by its wave function, a point tp 
in L^(r2,C). When the system is submitted to an electric 
field (e.g., a laser), the time evolution of the wave function 
is given, under the dipolar approximation and neglecting 
decoherence, by the Schrodinger bilinear equation: 



(-A + V{x))^P{x, t) + u{t)W{x)^{x, t) (1) 



where A is the Laplace-Beltrami operator on fi, V and W are 
real potential accounting for the properties of the free system 
and the control field respectively, while the real function of 
the time u accounts for the intensity of the laser 

In view of applications (for instance in NMR), it is 
important to know whether and how it is possible to choose 
a suitable control u : [0, T] R in order to steer ([T} 
from a given initial state to a given target. This question 
has raised considerable interest in the community in the last 
decade. After the negative results of [1] and [2] excluding 
exact controllability on the natural domain of the operator 
— A + 1/ when W is bounded, the first, and at this day the 
only one, description of the attainable set for an example of 
bilinear quantum system was obtained by ([3], [4]). Further 
investigations of the approximate controllability of ([U were 
conducted using Lyapunov techniques ([5], [6], [7], [8], [9], 
[10]) and geometric techniques ([11], [12]). 

B. Various notions of energies 

Quantum control is a trans-disciplinary field where differ- 
ent communities use the same word "energy" with possibly 
different meaning. 

Mathematically, the energy of system ([T]i is any norm in 
(a subspace of) C) and the energy for the control 



u in any norm in the space of admissible controls. A 
recurrent issue when studying systems of the type of ([l) 
is to obtain a priori estimates of the energy of the system in 
terms of some energy of the control. Such energy estimates 
are crucial for many reasons, both for mathematical and 
engineering purposes, including for instance the proof of 
the well-posedness of the system and the regularity of the 
solutions [13], or estimates of the distance between the 
original infinite dimensional systems and some its finite 
dimensional approximations (see Section Hl-CI below). 

Physically, the energy of the quantum system ^ with 
wave function is £'(V') — Jq [(—A + Vyip\ ip d/i. The 
physical energy is therefore constant in time whenever the 
control u is zero. When the control u is nonzero, and 
provided suitable regularity hypotheses, the energy evolves 



dE 
dt 



= 2u{t)Q (^j^ [(A + V)W\ Wi) Ai^ . (2) 



Note that the time derivative of the energy E at time t 
depends on the value u{t) of the intensity of the external 
field and on the wave function 'il^{t). 

A natural question is to relate the mathematical energy of 
the control with the physical energy of the system. 

Standard candidates for these estimates are the norms 

I|w||lp(o.t) = i^^Q 1 for some suitable p > 0. 

Indeed, many previous works addressed the problem of the 
optimal control of the system ([T]i for costs involving the 
norm of the control (see for instance [14] or [15]). The main 
reason for choosing the norm is the fact that the natural 
Hilbert structure of allows the use of the powerful tools 
of Hilbert optimization. It is common belief that there is a 
natural relation of the norm of u and the energy of the 
system. The note [16] showed that, in general, the L^-norm 
provides more information on the evolution of the system 
than other L^-norms for p > 1. 

C. Framework and notations 

To take advantage of the powerful tools of the theory of 
linear operators, we reformulate the bilinear dynamics dl) in 
more abstract framework. In the separable Hilbert space H, 



we consider the bilinear system 

^{t)^Am+uit)Bm (3) 

at 

where the (time independent) Hnear operators A and B 
satisfy some regularity assumptions. 
Assumption 1: The triple {A, B, ^) is such that 

1) A is skew-adjoint, possibly unbounded, on its domain 
DiA); 

2) —iA is positive; 

3) B is bounded relatively to A: there exist a and 6 in R 
such that \\Bip\\ < a\\Ai(;\\ + b\\ip\\; 

4) $ = is a Hilbert basis of H made of 
eigenvectors of A: for every j in N, there exists Xj 
in R such that A(pj = —iXj(f)j; 

If A and B satisfy Assumption 1 1131 we denote 

= inf{ae R | 36 e R 
for which \\B^\\ < a\\A^\\ + b^jW, Wijj G D{A)}. 

It is known that if {A, B, $) satisfies Assumption [T] then 
for every u : [0, T„] {-1/\\B\\a, 1/\\B\\a) with bounded 
variation, there exists a continuous mapping i i~> T" taking 
value in the unitary group U(iJ) of H such that, for every ip 
in D{A), t T^ip is differentiable almost everywhere and 
satsifies ^ for almost every t in (0, T]. For a proof of this 
well-posedness result, see [17] for a general theory of time 
dependent (non-necessarily skew-adjoint) Hamiltonians or 
[18] for an elementary proof adapted to the bilinear structure 
of©. 

Definition 1: Let {A, B, $) satisfy Assumption [T] The 
system {A, B) is approximately controllable if, for every 
tpo,ipi in the unit Hilbert sphere, for every £ > 0, there 
exists : [0, T^] R such that HT^'-^q - ^"111 < £■ 
The following sufficient criterion for approximate controlla- 
bility is the central result of [12], centered on the notion of 
non-degenerate (or non-resonant) transitions. 

Definition 2: Let {A,B,^) satisfy Assumption [T] A pair 
(j, k) of integers is a non-degenerate transition of {A, B, $) 
if (i) {(j)j,Bcj)k) ^ and (ii) for every {l,m) in N^, |Aj - 
Xk\ = |A; - A„i| implies (j,k) = {l,m) or {(j3i,B(j)„,) = 
or {j,k}f^{l,m} = 0. 

Definition 3: Let {A,B,^) satisfy Assumption [T] A sub- 
set S of is a non-degenerate chain of connectedness 
of [A, B, $) if (i) for every (j, k) in S, {j, k) is a non- 
degenerate transition of [A^B) and (ii) for every rci,rb in 
N, there exists a finite sequence Va ~ ro,ri, . . . ,rp — ri, in 
N such that, for every j < p — 1, {rj, rj+i) belongs to S. 

Proposition 1: Let {A, B, $) satisfy Assumption [T] If 
{A, B) admits a non-degenerate chain of connectedness, then 
(A, B) is approximately controllable. 

D. Main result 

The contribution of this note is to show the good agree- 
ment between the total variation of the control and the 
variation of the A-norm of the wave function. The A-norm, 
equal to \\A%Ij\\ for every ij) in D{A) is not equal, in general, 
to the energy However, if 4>j is an eigenvector 



of A with associated eigenvalue — iAj, then = |Aj| = 

We have to distinguish between the cases where B is 
bounded and when it is not. 

1 ) Bounded case: When B is bounded, the growth of the 
A norm of Tfip is at most linear with respect to the total 
variation of the control (see Section III-B.ll i. We present, in 
Section IIV-AI an example for which the growth is indeed 
linear. More precisely, we will show the following. 

Proposition 2: There exists {A, B, $) satisfying Assump- 
tion [T] with B bounded such that, for every M in R, there 
exists : [0, Tm] R with bounded variation with 
||AT^-0i|| > M and M > mTV[o,T,,]{^M). 

2) Unbounded case: When i? is unbounded, the growth of 
the A norm of T'^ip is at most exponential with respect to the 
total variation of the control (see Section |ri-B.2t . We present, 
in Section IIV-BI an example for which the growth is indeed 
exponential. More precisely, we will show the following. 

Proposition 3: There exists a triple {A, B, $) satisfying 
Assumption [T] with B unbounded such that, for every 
M large enough in R, there exists um ■ [0, Ti\i] — >■ 
R with bounded variation with || AT^JJJ(/)i|| > M and 

M > 4exp |^y|||S|UTV[o,TM](uM)j - 6. 

E. Content of the paper 

In Section HIl we review some classical estimates for the 
growth of I v4 1*" -norms of the wave function in terms of 
norms (Section Ill-Al l and total variation (Section III-Bb of 
the control. Some examples of use of these estimates for 
the approximation of the infinite dimensional system I© by 
its finite dimensional approximations are given in Section 
III-CI Section Hill is a quick survey of basic facts about 
averaging theory for finite dimensional bilinear systems. 
These convergence results will be instrumental in Section 
ITVl to prove Proposition |2] (Section lIV-Ab and Proposition |3] 
(Section HV-Bl) . 

II. SOME ENERGY ESTIMATES 

A. Weakness of estimates 

Let (A, _B, $) satisfy Assumption [T| and admit a non- 
degenerate chain of connectedness. For every r > 0, for 
every j, fc in N and e > we define k) as the set of 

functions u : [0,T„] ^ R in Li([0,T„]) n L''([0,r„]) such 
that ||Ty^(/)j — (j)k\\ < £■ We consider the quantity 

Cr{(t>j,4'k) = sup ( inf ||w||l'-(o,t„) I • 

This quantity is the infimum of the i''-norm of a control 
achieving approximate controllability. It clearly satisfies the 
triangle inequality. Next proposition states that Cr is a 
distance on the space of eigenlevels only when 7 =1. Its 
proof is given in [16]. 

Proposition 4: Ci is a distance on the set {4>j,j G N}. 
For r > 1, Cr is equal to zero on the set {(pj,j G N}. 
Proposition |4] illustrates various flaws of estimates of 
system Q. First, and contrary to the immediate intuition. 



LP norms with p > 1 (and in particular the norm) 
do not permit to distinguish among the energy levels of 

A. Precisely, if [A, B, $) admits a non-degenerate chain 
of connectedness, for every non empty open set V in Sh, 
there exists u : [0,r] (-1/a, 1/a) in LP([0,r]) with 
II'^I|lp([o,t]) arbitrarly small such that T^0i belongs to V, 
see [19]. 

While the norm allows to distinguish among the energy 
levels of A, the distance Ci depends only on B and the 
non-degenerate chains of connectedness of {A,B,^) (and 
not on the eigenvalues of A). For instance, the computation 
of Ci(0i, ^2) done in Section IV of [16] remains valid, 
with unchanged result, if one replaces A by ±i|A|'^ for any 
positive integer k. 

B. Estimates based on total variation 

The following estimates can be deduced from the general 
theory due to Kato [17]. They are valid in context much 
broader than Assumption [T] In particular, there is no need 
for H to admit a Hilbert basis made of eigenvectors of A. 

In the following we will impose u(0) — 0. This is always 
the case if one replaces A hy A + u{0)B. Moreover if 
{A, B, $) satisfies Assumption[T]then there exists and b E 
R such that {A + u{0)B — ib, B, $') satisfies Assumption [T] 
as well. 

1 ) Bounded case: 

Proposition 5: Let {A, B, $) satisfy Assumption [T] with 
B bounded. Then, for every u : [0, T] ^ R with bounded 
variation and u(0) = 0, for every j in N, ||AT^(/)j|| — 
\\A<j,,\\<2\\B\\TV[o,T]iu). 

Proof: Notice that if u(0) = then \u{T)\ < 
TV[Q T]{u)- Hence it is enough to prove for every j in N 
that |'|(A + u(r)S)T^0,|| - \\A^,\\ < ||i3||TV[o,T]M. 

Any bounded variation function can be approximated 
pointwise by a sequence of piecewise constant functions 
(wn)neN such that |u„| < |m| and TV[o,t](u„) < 
TVio,T]{u)- 

Following [18], we have that T^"0j T^0j. Thus it 
is sufficient to prove the statement for piecewise constant 
controls. The proof for piecewise constant controls follows 
from the estimate \\{A + uB) exp (t{A + uB))(j)\\ - \\A(j)\\ < 
Indeed for a piecwise constant function the associ- 
ated T" is a product of expt.i{A + UiB) for different values 
of Ui and t, . The details of the proof are similar to those of 
[18, Section 2]. ■ 

2) Unbounded case: 

Proposition 6: Let {A, B, $) satisfy Assumption [T] with 
B unbounded. Then, for every < S < 1, for every u : 
[0,r] ^ (-(1 - 5)/||B|U,(l - S)/\\B\\a) with bounded 
variation and u(0) = 0, for every tp in D{A), \\AT}^^\\ < 

el|B|UTV[o,Tl(u)/<5||^^||, 

The proof in the unbounded case, which can be found in 
[18, Proposition 3], follows the lines of the bounded case. 

C. Good Galerkin Approximations 

Assumption 2: The quadruple {A, B, $, k) is such that 
1) (A, B, $) satisfies Assumption [U 



2) —iA is positive; 

3) fc is a positive real number; 

4) for every u in R, the domains D{\A + uB\'') of \A + 
uBl*" and D{\A\'') of \A\'' coincide; 

5) there exists d,r in R, r < k such that Hi??/' 1 1 < 

|j4|'"-!/'|| for every ip; 

6) the supremum Ck{A,B) of the subset of R 
{|5R(|A|V,i3V'>l/l(l^lV»U e DQAl'^)} is fi- 
nite. 

For every N in N, we define the orthogonal projection 

j<N 

Definition 4: Let N € N. The Galerkin approximation of 
© of order N is the system in H 

i=(AW+?ii3W)x (En) 

where A^-^^^ — ttnAttn and _B(^) — ttnBttn are the 
compressions of A and B (respectively). 

We denote by X^^ n)(^' '■^^ propagator of ([S^. 

Definition 5: The system {A, B, $) admits a sequence of 
Good Galerkin Approximations (GGA in short), in time T € 
(0, +00], in a subspace D (with norm || • of H, in terms 
of a functional norm N{-) on a functional space U if, for 
any K,e > 0, for any ip in D, there exists in N such 
that, for any u in U, N(u) < K implies \\{X]^^ r^^{t,0) - 
TIq)'iP\\d < e for any t <T. 

Proposition 7: Let {A, B, $, fc) satisfy Assumption |2] 
Then {A, B, $) admits a sequence of good Galerkin approx- 
imations in infinite time, in D{A) in terms of norm for 
locally integrable controls. 

Last proposition is proved in [20] for piecewise constant 
controls. The generalization to controls follows from [18]. 

Propositions: Let d > 0, r < 1 and (A, 5,$) sat- 
isfy Assumption [U with ||i3V'll 5; '^11 l^rV-"!! for every ip 
in _D(|A|'"). Then {A,B,^) admits a sequence of good 
Galerkin approximations in infinite time, in D{A) in terms 
of TV + norm for controls with bounded variation. 

This proposition is proved in [18]. Notice, that if imposing 
u{0) =0 for the control term then the norm of the control 
over any finite time interval is bounded by a multiple of the 
total variation TV. 

III. PERIODIC CONTROL LAWS OF BILINEAR 
QUANTUM SYSTEMS 

A. Averaging theory 

The mathematical concept of averaging of dynamical 
systems was introduced more than a century ago and has 
now developed into a well-established theory, see for instance 
the books of Guckenheimer & Holmes [21], Bullo & Lewis 
[22] or Sanders, Verhulst & Murdock [23]. It was observed 
that, for regular F and small e, the trajectories of the system 
X — eF{x,t,£) remain e close, for time of order 1/e, to 
the trajectories of the average system x ~ F{x) where 
F{x)=Viuit^^l/t J^F{x,t,0). 



In quantum physics, this concept of averaging is used 
intensively to transfer a system of type (|3) from an eigenstate 
of A associated with eigenvalue — iAj to another associated 
with eigenvalue — iAfc with a periodic control with small 
enough amplitude and frequency \Xj — Afc|. 

The following results is proved in [24]. 

Proposition 9: Let {A, B, $) satisfy Assumption [T] As- 
sume that {j, k) is a non-degenerate transition of {A, B, $). 
Define T = 2tt/\Xj - Xk\ and let u* : R ^ 

(-I/II-BIU, I/II-BIU) be T-periodic and with bounded 

rT 

variation on [0,r]. If / u*(T)e'(^^~^'='"'dr 7^ and 



/ M*(r)e'(^'"^'"'^dr = for every (Lm) such that (i) 

{Ik} ^ {l,m}, (ii) {j, fc} n {l,m} ^ 0, (iii) \Xi-\,n\e 
(N \ {l})|Aj — Afel and (iv) bi„i ^ 0, then, for every n 
in N, there exists T* in {nT* - T, nT* + T) such that 
I {(t>k , jY) (-^n ' I tends to 1 as 71 tends to infinity, with 



T* = 



216 



K 



IT* 



\u*iT)\dT, 



and C 



sup 

(j,fe)eA 



Jo " 



(T)e'(^'--^'")^dT 



sm vr 



■ |A,-A„j 



where A is the set of all pairs {I, m) in {1, . . . , A^}^ such that 
^ Oand {l,m}r^{i,k} ^ and |A,-A„,| i Z|A2-Ai|. 
Notice that Proposition |9] does not claim that 
||A(A'"^ ^j(T*, 0)(/ij — 4>k)\\ tends to zero as n tends 
to infinity. However, 

liminfpX(V^)(r,:,0)(/.,|| 

> liminf Afc|((/.fe,X(\^)(r„*,O)0,)| = A^. 

Using Propostion|7]or Proposition |8]these can be extended 
to inifinite dimensional system ^ with {A, B, $, fc) satisfing 
Assumption I2] or Assumption [U with !i '^Ill^rV'll for 

every ip in _D(|^|'') for some d > and r < 1. 

B. Averaging using the sine function 

Let (A, B, $) satisfy Assumption [T] and (j, fc) be a non- 
degenerate transition of (A, B, $). We define w = |Aj — 
Xk\- We apply Proposition |9] with it* : i i-^- sin(wt). For 
71 large enough, \\u*{t)/n\\Loa < Straightforward 
computations give 



r = — , T* 

UJ 



I 



and we compute 




uj \ cos{ujt)\dt 



co8{ujt)\dt' 



I cos{u}t)\dt 



As 77 tends to infinity, — 

77, 

hence 



I cos{ujt)\dt tends to zero, 



lim TV[o,T*i I — 



n 



lim — 

00 77 







T 


L 



I coB{ujt)\dt 



lim 4 



2uj 

\bjk\ 



IV. EXAMPLES 
A. The bounded case: 2D rotation of a linear molecule 

Consider a linear molecule whose only degree of freedom 
is the planar rotation, in a fixed plan, about its fixed center 
of mass. This system has been thoroughly studied (see the 
references given in [25] or [20] for instance). 

In this model, the Schrodinger equation reads 



•^^ A / 



(4) 



il ~ R/27rZ is the unit circle endowed with the Riemannian 
structure inherited from R, H is the space of odd functions of 
L^(ri, C), A = iA (A is the restriction to H of the Laplace- 
Beltrami operator of fl) and B : ^ [9 cos{9)tp{9)) is 
the multiplication by cosine. 

In the Hilbert basis ^ = {9 sm{k9))ke-N of H, A 
is diagonal with diagonal — ifc^, k = 1 ... 00 and B is tri- 
diagonal with bk,k — 0,bk.k+i = bj_k = for every 

k,j in N such that \j — A;| > 1. 

The triple {A, B, <&) satisfies Assumption[Tl B is bounded 
and \\Bijj\\ < 0\\Aip\\ + V2\\^p\\ for every ^ in H. The 
set {{k,k + l),fc e N} is a non-degenerate chain of 
connectedness for {A, B, $). 

For every j, nj in N, we define the control u*'^'"^ : t G 



[0, 2r7j7r] I-)- sin((2j 



l)t)/77j, and for every N in N, 



we define -u*:("i."2:. .,niv-i) \jy (he concatenation of u 



i*,2,n2 



By Proposition |9] 

liminf ||AT"*'*"''"' 

From Section IIII-B[ we compute 

liminf TVro,27r(ni+n2 + ...+n;N 
Til ,712 ,. . . ,njv — 1— >-oo 



> A 



N 



*,(T^l,n2,...,nJv-l)^ 



N-1 

4(2j + 1) = 47V2 



which proves Proposition |2] 



B. The unbounded case: perturbation of the harmonic oscil- 
lator 

The second model we consider is a perturbation of the 
quantum harmonic oscillator, with dynamics give by 



dt 



[(-A + x^) + (-A + a;2)-i] ^ + u{t)x^ip. (5) 



With the notations of Section |LC] H is the Hilbert 
space of the odd functions of L^(R, C), A = 
— i [(—A + x^) + (—A + x^)"^] where A is the restriction 



of the Laplacian to the space of odd functions and B is 
the multiplication, in H by 



~ix . Denoting by iJ„ the 
Hermite function, we check that AH2n-i = — i((4n — 
1) + (4n - iy^)H2n-i, hence $ = (i?2n-i)neN is a 
Hilbert basis of H made of eigenvectors of A. Moreover, 
BHi = -i{l/2Hi + ^/3/2H3) and, for every n in N, n > 2, 



BH2n-l — — i 




In the basis $, A is diagonal with diagonal entries (— i((4n— 
1) + (4n — l)~^))„gN and B is tri-diagonal. The system 
{A,B,^) is tri-diagonal in the sense of [18] and satisfies 
Assumption [U with \\B\\a < ^ (Proposition 12 of [18] 
applied with r — 1 and C — 1/4). 

For every j,nj in N, we define the control u*'^'"^' : t G 
[0,2nj7r] n- sin((2j + l)t)/nj, and for every N in N, 
we define y*.("i,ri2,.. i^y tjjg concatenation of 

y*, 2,712 y*,-'V-l.nN-l 

By Proposition |9] 

liminf HAT"*''"""' """'Vill > Ajv > 4Ar- 1. 

ni,...,n7v — 1 — ^oo 



From Section IIII-BI we compute, similarly to what we 
have done in Section Hl-B. II 

liminf TVio.2.(„,+„,+...+„„_,]K'("^'"^'-'""-^^) 

ni, n2, — 1~^CX3 

N-1 

= T- 

N 

< 21og(7V + l), 

which proves Proposition |3] 



V. CONCLUSIONS 



A. Contribution 



We exhibited two examples showing that the Kato esti- 
mates for the A-norm of the solutions of a bilinear quantum 
system, with bounded or unbounded coupling term, are 
optimal up to a multiplicative constant. These estimates are 
given in terms of the total variation of the control. 

B. Perspectives 

An interesting, and probably difficult, question is the 
optimal control of bilinear quantum systems when the cost is 
the total variation of the control. Approximation procedures 
(as the Good Galerkin Approximations presented in this note) 
allow to consider only a finite dimensional problem. The 
main difficulty will come from the non-smoothness of the 
cost (total variation) which will lead to the use of tools of 
non-smooth analysis. 
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